Fluctuating parts of nuclear ground state correlation energies 
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Background: Heavy atomic nuclei are often described using the Hartree-Fock-Bogoliubov (HFB) 
method. In principle, this approach takes into account Pauli effects and pairing correlations while 
other correlation effects are mimicked through the use of effective density-dependent interactions. 

Purpose: Investigate the influence of higher order correlation effects on nuclear binding energies 
using Skyrme's effective interaction. 

Methods: A cut-off in relative momenta is introduced in order to remove ultraviolet divergences 
caused by the zero-range character of the interaction. Corrections to binding energies are then 
calculated using the quasiparticle-random-phase approximation (QRPA) and second order many- 
body perturbation theory (MBPT2). 

Result: Contributions to the correlation energies are evaluated for several isotopic chains and 
an attempt is made to disentangle which parts give rise to fluctuations that may be difficult to 
incorporate on the HFB level. The dependence of the results on the cut-off is also investigated. 

Conclusions: The improved interaction allows explicit summations of perturbation series which 
is useful for the description of some nuclear observables. However, refits of the interaction parameters 
are needed to obtain more quantitative results. 

PACS numbers: 21.10.Dr, 21.30.Fe, 21.60.Jz 
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I. INTRODUCTION 

The atomic nucleus is a complicated quantum mechan- 
ical system where the probability to find a nucleon in a 
certain position is a function of the positions of the other 
nucleons. This is generally referred to as the nucleons 
being correlated and makes the wave functions of heavy 
nuclei too complex to compute directly using ab initio 
theory. One therefore has to resort to more tractable 
methods which take the most important correlation ef- 
fects explicitly into account, i.e. the ones that are impor- 
tant in order to describe observables, while the remaining 
effects are treated in an approximate way. 

An often used starting point is to assume that the in- 
medium interaction between nucleons can be modeled us- 
ing effective density-dependent internucleon potentials. 
Such potentials are generally employed in Hartree-Fock- 
Bogoliubov (HFB) calculations and their parameters are 
fitted to reproduce a number of experimentally known 
data on individual nuclei and to what is known about 
nuclear matter. In this way, Pauli effects and pairing 
type correlations are taken into account explicitly while 
the effects of other of correlations are described in an 
average way. This approach has the great advantage of 
being applicable to the entire nuclear chart at a reason- 
able computational cost. In the quest of more accurate 
nuclear models an important task however, is the sys- 
tematic investigation of which type of correlation effects 
can be modeled with the HFB method and which need a 
more explicit treatment. 
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Several studies have shown that going beyond a HFB 
treatment and adding corrections to binding energies re- 
sulting from shape vibrations, especially of quadrupole 
type, give an improved description of experiment [l|-|j|- 
These corrections are often taken into account by ei- 
ther using the generator-coordinate method (GCM) [l[ 
or through the use of a collective model e.g. a Bohr 
Hamiltonian Q. 

Alternatively, many-body perturbation theory 
(MBPT) offers a way to explicitly and pictorially include 
elementary processes that one might suspect to be 
responsible for correlations in different systems. For 
instance, within the so called random-phase approx- 
imation (RPA) one allows for an infinite number of 
particle-hole pairs to be excited out of the Hartree-Fock 
ground state and for multiple scattering between excited 
particles and holes. If the excitations instead are made 
out of the HFB ground state the same approximation is 
referred to as the QRPA designating RPA for Bogoliubov 
quasi-particles. An even simpler step beyond HFB is the 
second-order many-body perturbation theory (MBPT2) 
starting from the HFB ground state. Clearly, the virtual 
excitations included in this approximation is a subset of 
those included in the QRPA and in the present work we 
show results from both levels of approximation. 

Most of the effective nucleon potentials involve contact 
terms, i.e, interactions of zero range. This is certainly the 
case for interactions of the Skyrme type and such inter- 
actions give rise to divergences when going beyond the 
HF level. This can be seen e.g. by solving the two-body 
problem for 2 H analytically using contact interactions. 
Then the resulting binding energies become infinite [H, Q . 
The two methods used in this work for going beyond the 
HFB level include infinite summations of intermediate 
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states which inevitably leads to the same divergences in 
connection with zero-range forces. 

It is a major theme of the present work to eliminate 
such divergences by introducing cut-offs in momenta for 
our chosen Skyrme-like interaction potentials 0|. This 
procedure implicitly assumes that structures in binding 
energies as functions of nucleon number originate in cor- 
relation effects caused by the low-momentum part of the 
internucleon forces. And according to the results of the 
present work this assumption does not appear to be that 
far fetched. 

This paper is organized as follows: in Sec. II the reg- 
ularized Skyrme interaction is introduced. In Sec III 
we discuss the treatment of correlation effects using the 
Quasiparticle-random-phase approximation (QRPA) and 
the MBPT2 method. In sects. IV and V we analyze and 
discuss the results of our calculations. 



II. LOW-MOMENTUM INTERACTION 

A. Two-body interaction in the particle-hole 
channel 

A general two-body interaction that preserves the cen- 
ter of mass coordinate of the interacting particles can be 
expressed as 

V (r' ir ' 2 rir 2 ) = v (r',r)S(R — R') , 

where R = ^ (ri + r 2 ) and r = r\ —r 2 denotes the center 
of mass and relative coordinates respectively. The part 
of the potential depending on relative coordinates can 
be transformed to momentum space and for this part we 
adopt Skyrme's expansion Q given by 



v (k' k) = 7T 
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+ t 2 (1 + x 2 P a ) k' k + iW a (<ri + a 2 ) k' xk]. 

In this expression we have omitted the tensor potential 
included in Ref. @ since it is not used in the parametriza- 
tions we will employ later. This expression can be viewed 
as the first terms in a low-momentum expansion of the 
effective nuclear potential going up to second order in 
relative momenta (compare e.g. [9Mll| for higher or- 
der expansions). Although this form gives a reasonable 
description of the low-momentum parts, the expansion 
becomes unrealistic for large momentum transfers. As 
we will demonstrate later, for Hartree-Fock calculations, 
only the low-momentum matrix elements are important 
and the unphysical contributions generated by the ex- 
pansion for higher momenta can be ignored. 

For studies beyond the mean-field level however, the 
interaction gives diverging results unless some kind of 



truncation is enforced e.g. a truncation in excitation en- 
ergy. Nevertheless, in some beyond mean field calcula- 
tions, such as QRPA calculations, the results for low- 
lying states [IH and giant resonances [l3[ are in reason- 
able agreement with experiment indicating that the in- 
teraction may indeed have a wider applicability beyond 
purely mean- field calculations 1 . 

In order to investigate how well higher order correc- 
tions can be described using the low-momentum part of 
Skyrme's interaction we follow Skyrme's original sugges- 
tion [f| and introduce a cut-off in momenta. We replace 
his original interaction by 



,(A) 



(k',k) = v(k',k)9(A-k')9(A-k), (I) 



which vanishes at momenta above A (fm) _1 . In the limit 
of a large A one regains the results of the original untrun- 
cated interaction, but for finite values, the cut-off regu- 
larizes the interaction so that beyond-mean field calcula- 
tions converge. The introduction of the cut-off destroys 
the nice analytical properties of the zero-range interac- 
tion and increases the computational cost of calculating 
matrix elements. 



B. Two-body interaction in the particle- particle 
channel 



In the pairing channel we use the same finite-range 
separable-Gaussian interaction as was used in our previ- 
ous studies (l2l [l3j . Since this interaction has a finite 
range, no regularization is needed. We adopt an isospin 
invariant form, active in the T = 1 channel and use the 
same range parameter (a = 0.66 fm) as before. Since 
we will only consider cases where neutrons are in open 
shells we tune the pairing strengths to make the lowest 
neutron quasi-particle energies to agree with the exper- 
imental gaps determined in The resulting isovec- 
tor pairing strength becomes 560 MeV(fm)" 1 when the 
SKX Skyrme parameters [l5[ are used in the particle-hole 
channel and somewhat larger (640 MeV(fm) -1 ) when the 
SLy5 parameters [l6| are used. 



Density-dependent part of the particle-hole 
interaction 



Skyrme's expansion of the two-body potential is often 
used together with a density-dependent zero-range poten- 
tial which is intended to describe missing three-body and 



1 Note that in the calculations of low-energy excitations the dis- 
cussed divergences did not constitute a problem. In fact, it is 
mainly the high-energy excitations which are modified by a mo- 
mentum cut-off. 
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higher order contributions as well as giving a simple rep- 
resentation of missing many-body effects. The density- 
dependent terms cause difficulties when going beyond the 
mean field and different recipes to define a residual in- 
teraction exist in the literature [l7| . In this work we 
are mainly motivated by the success of the Skyrme in- 
teraction in connection with RPA type calculations and 
hence define the residual interaction as the so-called RPA 
residual interaction using the second derivative of the HF 
energy [l8| 
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The use of the RPA residual interaction in configura- 
tion interaction type calculations has been thoroughly 
discussed and investigated before using a non-regularized 
Skyrme interaction [19j. 

The density-dependent two-body interaction is intro- 
duced in the standard form 201 



vx-' (r',r)(p(R)) a S(R-R'), 



with a dependence on the nucleon density p to some 
power a which take on different values for different pa- 
rameterizations. The part dependent on relative coordi- 
nates is expanded to lowest order in relative momenta 
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and regularized with the same cut-off procedure (A- 
truncation) as used for the density-independent parts. 

In the practical calculations of matrix elements we 
start from a spherical Harmonic-oscillator basis and 
transform the basis functions to momentum space. The 
A-truncation can then be implemented using the Moshin- 
sky transformation [2f| to transform the coupled two- 
particle states to functions of relative and total mo- 
menta. Finally we employ the Pandya transformation 
f22l | to obtain matrix elements in the particle-hole channel 
and use the Wigner-Eckart theorem to obtain angular- 
momentum reduced expressions. The full implementa- 
tion of this new regularized potential was done by ex- 
tending the program HOSPHE (vl.02) (23| . 

One of the interactions employed in this work (SLy5) 
uses the direct part of the Coulomb interaction together 
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Figure 1: (Color online) The energy per nucleon in symmet- 
ric nuclear matter (equal numbers of neutrons and protons) 
shown for the SLy5 |la | interaction and for different A (fm) -1 
values in the Hartree-Fock approximation. The full line il- 
lustrates the standard result for the untruncated interaction 
(A = 00). 



with a Slater approximation for the Coulomb exchange. 
The Slater approximation results in a density-dependent 
term which mimics the HF exchange energy. In order to 
treat the HF part and the additional residual interaction 
consistently we have regularized the Slater term in the 
same way as for the other parts of the interaction. 

In order to have a first idea about the influence of 
the A -truncation we consider isospin-symmetric nuclear 
matter in the Hartree-Fock approximation. The corre- 
sponding zero-temperature equation of state (EOS) i.e. 
the energy per nucleon as a function of density is shown 
in Fig. [1] As seen in this figure, a A value of «1.5 (fm) -1 
leaves the EOS unchanged up to about twice the satura- 
tion density and a value of sal. 75 (fm) -1 leaves the EOS 
unchanged up to about three times the saturation den- 
sity. Thus, when choosing a value for the regularization 
we will consider values above f .5 (fm) -1 which keeps the 
relevant part of the EOS approximately the same. 



III. RESULTS FOR THE TOTAL 
CORRELATION CONTRIBUTION 

An expression for the RPA correlation energy in the 
quasi-boson approximation (QBA) was derived in (l8| 
and using an analogous derivation one obtains a corre- 
sponding expression in the QRPA case [24| 



Eqrpa = - ^2 ^ \ Y w 1 



(2) 



k<k' 



In order to evaluate this expression we start by defining 
matrices containing the positive energy QRPA vectors 

X= [X\X 2 ,...,X N ] ,Y= [Y\Y 2 ,...,Y N ] 
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and corresponding energies 
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Then starting from the QRPA equation [22 



A B \ (X 
B* A* [ Y 



one can write the pair of equations 



X 

-Y 



YX^A + YX^BYX- 1 = YnX- 1 

B*+A*YX~ 1 = -YflX- 1 . 

Summing these equations together and introducing C 
YX" 1 the result is the equation 



B* + A*C + CA + CBC = 



(4) 



which is similar to the multiple scattering series derived 
in Ref. (2a |. In terms of C, the QRPA correlation energy 
becomes 
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Figure 2: (Color online) The upper panel shows the binding 
energy contribution from QRPA and MBPT2 scaled with the 
HFB energy. The lower panel shows the influence of the trun- 
cation on the HFB energy. Calculations were done using the 
SLy5 [T?3 | parametrization for the Skyrme interaction using a 
basis of 14 oscillator shells. 



Eqrpa — - ^2 ^kk' .ll'Ckk' ,ll< ■ 
k<k' ,KV 

Furthermore, by splitting the A matrix 

Akk',w = (Ek + Ek>) SkiSk'i' + Akk',w 
Eq. [4] can be written 

-1 



(5) 



Ckk'W — 



( B k',w + (A*C) kk , w + (CA) kk , M , 
(CBC) kk , a ,) . (6) 



Finally, assuming that an iteration procedure for C 
converges, we can evaluate C order by order where the 
first order contribution 
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is obtained by putting C equal to zero on the right hand 
side of Eq. [6] Higher orders C^ n > are thus obtained by 
repeatedly inserting the previous expression G 1 ( Tl-1 - ) on 
the right hand side. By using in the formula for the 
correlation energy, Eq. [5j we obtain Eg'j i p A . Numeri- 
cally we have also verified that the iteration converges to 
results consistent with Eq. [2] 

An alternative approach which is less costly numeri- 
cally is to evaluate the correlation contribution from sec- 
ond order perturbation theory [2(| starting from the HFB 



ground state and treating the residual part of the quasi- 
particle Hamiltonian as a perturbation. In this case, the 
only contributions that arise come from scattering to four 
quasi-particle states via the H 4a part [l8| of the Hamil- 
tonian. This contribution can be expressed in terms of 
the QRPA B matrix according to 



E 



(2) 

MBPT 



= 1 E 



\B 



ll'kk' 



6 ^ E k + E k > 
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E v + Ev 
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It is interesting to compare the QRPA and MBPT se- 
ries order by order. The lowest order QRPA term is three 



times larger than E 



(2) 

MBPT 



while the third order E 



(3) 

MBPT 



is exactly obtained in the QRPA series. In higher or- 
ders, the two series differ as the QRPA expression only 
includes a subsequence of the full MBPT series. 

In earlier papers by Ellis [13, [H|, the RPA correla- 
tion energy was investigated by starting from an un- 
paired ground state and summing contributions from 
both normal and pairing vibrations using diagrammatic 
techniques. In this way it was shown that in the QBA, 
the second order contribution appears twice in the sum- 
mation of the particle-hole ring series and once in the 
particle-particle series. A suggested remedy for this 
overcounting was to remove the second-order term from 
the particle-hole series and only keep it in the particle- 
particle series. In this work we have however refrained 
from using this approach since it is not directly applica- 
ble when starting from a HFB state where normal and 
pairing vibrations are generally mixed. 

The corresponding correlation energies evaluated with 
the two methods described above are shown in panel (a) 
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of Fig. [21 In both results we have not included the part 
of the B matrix obtained in pnQRPA [22| which is associ- 
ated with excitations of proton- neutron pairs. Although 
this contribution is certainly interesting, a first step in the 
direction of including these effects would involve tuning 
the effective interactions in the T = pairing channel. 

As seen in Fig. [21 for both methods the correlation 
energy amounts to a rather large part of the total binding 
energy. The QRPA formula predicts the largest values 
as the QBA overestimates the ground state correlations 
(l8l , [2tN-29| | . Although this could possibly be corrected 
for, in the following we will instead focus on the MBPT2 
results. 

In the case of 16 0, the smallest A = 1.4 (fm) _1 used 
in the figure gives a contribution from MBPT2 which is 
21 % of the HFB energy. This contribution gradually 
decreases for the heavier nuclei and becomes 16 % in 
132 Sn. 

Panel (b) of Fig. [2] shows the influence the regular- 
ization has on the HFB energy. As seen in this figure, 
the HFB energy converges to the untruncated value as 
the cut-off is increased and even for rather low cutoffs of 
A = 1.5 (fm) _1 the change in total binding energy stays 
within a few percent. This tells us that the HFB energy 
is not very sensitive to higher momentum parts of the 
potential in the particle-hole channel. 

In both the QRPA and in MBPT2, the correlation en- 
ergy increases rapidly with increasing A and in the fol- 
lowing we will consider A values in the range of 1.6-1.8 
(fin) -1 . These values lead to the smallest correlation 
energies while causing moderate changes of the HFB en- 
ergies. 

Since the HFB energy stays roughly the same, it is a 
good approximation to neglect the regularization for the 
HFB part of the calculation and only regularize when 
generating the residual interaction. This is quite im- 
portant in order for the method to be practical since 
otherwise one would have to recalculate the regular- 
ized density-dependent interaction in each HFB iteration. 
This is a strong motivation for introducing the regulariza- 
tion in the way done here rather than using for example 
a Gaussian interaction. For the purpose of making the 
least amount of approximations, we have however used 
the time-consuming strategy to recalculate the regular- 
ization in each HFB step. 

Using an angular-momentum coupled notation, the 
correlation energy can be divided into partial contribu- 
tions arising from QRPA excitations with different total 
angular momentum J and parity 7r. Since the MBPT2 
result can be seen as an approximation to the full QRPA 
results we use the same division into multipole contri- 
butions also in this case. These partial contribution 
are shown in Fig. [3] for 132 Sn. As seen in this figure 
the largest contributions come from natural parity states 
with (— l)' 7 = 7T. Both positive parity and negative par- 
ity contributions are equally important and show maxima 
for J — 4 and J = 5 respectively. 
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Figure 3: (Color online) Partial contributions to the MBPT2 
correlation energy for A = 1.8 and 132 Sn using SLy5. 



IV. FLUCTUATING PARTS OF THE 
CORRELATION ENERGY 

The internucleon potentials we employ have parame- 
ters which are fitted in order to give reasonable nuclear 
properties at the HFB level. Therefore the average part 
of the correlation energies is already effectively included 
through the fitting of the model parameters. The to- 
tal ground-state energy can be divided into a liquid-drop 
part that captures the average variations as a function 
of nucleon numbers and a fluctuating part that mainly 
depends on the shell structure. Rather than perform- 
ing a full refit of the interaction to have a model on the 
MBPT2 level, in this first study we will make a simple 
compensation for this overbinding. We compensate by 
fitting a liquid-drop expression |18| 

Eld = a vo iA + a sur fA 2 ^ 3 + a sym — (8) 

(2) 

to the correlation energy Ew BPT which is then sub- 
tracted to give the fluctuating part of the correlation en- 

(2) 

ergy AE = E MBPT — Eld- The main goal is to get an 
idea what kind of fluctuations one obtains and to see if 
these are correlated with errors obtained in the descrip- 
tion of ground-state energies. 

In this way the renormalized correlation contributions 
associated with 2 + and 3~ vibrations are extracted and 
shown in Figs. [4] and [5j The fluctuations show pro- 
nounced shell effects and tend to give increased bind- 
ing energy contributions for the open shell nuclei com- 
pared to the magic ones. The two different choices of 
A-truncation shown in the figures give similar results in- 
dicating that the obtained fluctuations are mainly asso- 
ciated with the properties of the low-momentum part of 
the interaction. 

In the 2 + channel, the SKX interaction gives larger 
fluctuations than the SLy5 interaction. Previous results 
for quadrupole correlations using the SkI3 interaction 
and a collective Hamiltonian f30j | gave similar results with 
fluctuations that are somewhere in between the ones we 
get for the SKX and the SLy5 interactions. 
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Figure 4: (Color online) The renormalized part of the MBPT2 
correlation energy associated with quadrupole shape vibra- 
tions (J = 2 + ) shown for A = 1.6 (dashed curve) and A = 1.8 
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Figure 6: (Color online) The renormalized part of the MBPT2 
correlation energy separated into contributions from different 
multipoles for A = 1.8. Solid curves are for SLy5 and dashed 
curves for SKX. In this figure the renormalization was done 
by fitting Eq. [8] to Sn nuclei only. 
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Figure 5: (Color online) Same as Fig. [4] for octupole shape 
vibrations (J — 3~). 
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Figure 7: (Color online) Difference between experimental 
[3ll . l32l | and theoretical ground-state energies in the HFB ap- 
proach (dotted curves) and when adding surface vibrations 
corresponding to multipolarities J n — + , 2 + , l - and 3~ (full 
curves). Figure is drawn using the SLy5 interaction and 
A = 1.8 (fm)" 1 . 



Contributions from octupole vibrations are similar in 
magnitude to the quadrupole vibrations and show the 
same tendency of increasing the energies for magic nuclei 
as compared to their neighbors. Notable exceptions are 
16 O and 40 Ca which have Fermi levels between opposite 
parity shells and show the reversed trend. 

Fig. [HI shows the fluctuating part of the correlation 
energies separated into contributions from different mul- 
tipoles in the case of the Sn chain. The main fluctuations 
come from the J T = 2 + , 4 + , 3~ and 5~ channels and show 
the trend of making double-magic nuclei less bound rel- 
ative to the semi-magic ones. Interesting exceptions are 
found in the 10 + and 11~ contributions which show the 
opposite trend around N=82. Going even higher in mul- 
tipoles the curves tend to flatten out. 



A. Comparison with experiment 

When comparing with experimental ground-state en- 
ergies, we restrict ourselves to contributions from the 
well studied low-order multipoles J x = + ,2 + ,l~ and 
3 _ where the effective interactions generally give rea- 
sonable results for low-lying collective states and giant 
resonances. 

The difference between calculated and experimental 
ground-state energies using the SLy5 and the SKX [l5| 
interactions are shown in Figs. [7] and [3] The dotted lines 
denote the results of the HFB treatment (A = oo) using 
a total of 31 oscillator shells. For both interactions, the 
HFB results are within a few MeV of the experimental 
values. The SKX results differ somewhat from the ones in 
[l5| which is due to a different treatment of the Coulomb 
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Figure 8: (Color online) Same as [7] but for the SKX interac- 
tion. 



interaction. While we calculate the direct Coulomb con- 
tribution directly from the proton density, in [l5T | some 
additional corrections were taken into account. 

It interesting to notice that with both interactions the 
errors for magic nuclei with TV = 8, 20 and 28 go up in 
energy compared to their neighbors which could possibly 
be cured by making the corresponding gaps in the neu- 
tron spectra somewhat larger. With SLy5, the situation 
is reversed for the gaps at TV = 50 and 82 where the 
errors instead dip down. 

The lowest order surface vibrations are calculated us- 
ing MBPT2 and a total of 14 oscillator shells. The fluc- 
tuating part of the correlation energy is well converged 
using 14 oscillator shells and is extracted by fitting the 
liquid-drop expression to all the included nuclei. The 
solid lines in Figs. [7] and [S] show the results of adding 
these fluctuating parts to the HFB energies. As discussed 
previously the main effect of the surface vibrations is to 
push the magic nuclei up in energy as compared to the 
neighboring nuclei. For SLy5 this gives corrections that 
go in the right direction in the region of the TV = 50 and 
82 gaps but in the opposite direction for the lighter magic 
nuclei. Although both interactions give rise to similar 
fluctuations, the larger magnitude fluctuations in com- 
bination with different HFB results obtained using the 
SKX parametrization, compares less favorable with ex- 
periment. If a Skyrme interaction tuned at the MBPT2 
level is used on the HFB level, one would expect it to 
predict the magic nuclei to be more bound than their 
neighbors in order to leave room for the additional cor- 
relation part. 

Some of the nuclei included in the plot have the same 
number of neutrons and protons which gives rise to an 
additional contribution to the binding energy. This con- 
tribution is often modeled by adding so called Wigner 
corrections (see e.g. [33[) which gives an additional bind- 
ing energy contribution of roughly 2 MeV for the N = Z 
nuclei. Such a contribution would reduce some of the 
remaining fluctuations but would not improve the re- 
sults around 48 Ca for example. Furthermore, such a phe- 
nomenological treatment is clearly unsatisfactory and a 



more thorough investigation of these interesting effects is 
clearly called for. 

In the case of A = 1.8 (fm)" 1 the parameters obtained 
from the liquid-drop fit to the SLy5 correlation energies 
resulting from J = + , 2 + , 1 _ and 3 _ vibrations become 
{a vo i,a surf ,a sym } = {0.99,-8.24,0.77} MeV, while for 
SKX the average contribution is roughly twice as large. 
A possible reason that the SKX interaction gives more 
correlation energy is that SKX has larger effective mass 
(m*/m = 0.99) [H than SLy5 (m*/m = 0.69) QJ] and 
thus a denser spectrum, giving smaller denominators in 
Eq.0 

Typical liquid-drop parameters obtained when 
fitting to experimental ground-state energies are 
{a vo i,a surf ,a sym } = {-15.68,18.56,28.1} MeV 
Thus, the main change in the average energy obtained 
by adding the correlations resulting from low-lying 
surface vibrations is to modify the surface energy. The 
reduction of the surface energy and increased energy for 
the volume part will likely move nucleons from the bulk 
to the surface leading to a more diffuse surface region. 
Thus refitting the Skyrme parameters to absorb the 
average part of the correlations and have a model on 
the MBPT2 level would likely involve timing not only 
the density-dependent terms but also the gradient terms 
which are more sensitive to the surface region. When 
refitting, it is important to have as small correlation 
corrections as possible, so that the HFB ground state 
is a reasonable first approximation. In this respect 
the smaller average contribution obtained in the SLy5 
case makes it a better starting point. Nuclear matter 
properties can also be used to refit, but then the nuclear 
matter EOS has to be calculated at the corresponding 
level of many-body theory (see e.g. [34| for a description 
of nuclear matter at MBPT2 order). 



V. SUMMARY AND CONCLUSIONS 

The problem we set out to investigate was whether 
effective nuclear interactions can provide improved de- 
scriptions of nuclear binding energies when correlation 
effects beyond the HFB level are taken into account. To 
this end, it was essential to introduce a momentum cut- 
off in Skyrme's potential in order to obtain convergent 
results. The calculations show that even with a low cut- 
off, the average part of the correlation corrections are 
quite substantial (about 25 % of the total binding energy 
with MBPT2 and A = 1.8 (fin) -1 ). We then considered a 
schematic renormalization by removing the average parts 
of the correlation energies. The remaining fluctuations 
are similar for both interactions studied and not so sen- 
sitive to the exact choice of the momentum truncation. 
When the SLy5 Skyrme parametrization is used, the fluc- 
tuations associated with low-lying surface vibrations does 
lead to a reduction of the errors compared to experiment. 
In order to obtain a model that can be used with more 
confidence a refit of the interaction parameters should 
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be performed. The ideal would be to compare results of 
an interaction fitted on HFB level to those of an inter- 
action fitted on the MBPT2 level using the same set of 
experimental data. 

Some interesting features can anyway be learned from 
the obtained fluctuations. One result is that octupole vi- 
brations are predicted to give fluctuations of similar mag- 
nitude as the quadrupole corrections and to contribute in 
a similar way. It is also interesting to see that higher mul- 
tipoles such as 4 + and 5~ gave rise to large fluctuations 
in the case of Sn isotopes. 

The fluctuating parts were extracted using MBPT2 but 
the QRPA formula is also promising in the sense that it 
allows an infinite summation of diagrams. However, in 
order for it to be a practical tool, the QBA approximation 
must be improved and a careful study of the corrections 
in the quasi-particle case would be needed. Once such 
a formalism is in place, the correlation energy could be 
calculated using iterative approaches (35| similar in spirit 
to the ones we recently employed for the calculation of 
low- lying excitations [12j ■ 

In summary, we have regularized Skyrme's potential 
and used it to study higher order corrections to binding 
energies beyond the HFB approach. Compared to other 



approaches, the method used here has the advantage of 
not relying on energy truncations in order to converge 
and that correlations resulting from many degrees of free- 
dom (e.g. vibrational modes) can be simultaneously in- 
cluded. Apart from nuclear binding energies studied in 
this work, there are other quantities that could possi- 
bly be modeled better in a formalism that goes beyond 
the HFB approximation. An example is the calculation 
of alpha-decay preformation amplitudes which shows a 
dramatic increase as correlations between nucleons are 
introduced (36| . 
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